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Abstract. We analyze the structural and dynamical properties of disk-like 
objects formed in fully consistent cosmological simulations which include 
inefficient star formation. Comparison with data of similar observable prop- 
erties of spiral galaxies gives satisfactory agreement, in contrast with pre- 
vious findings using other codes. This suggests that the stellar formation 
implementation used has allowed the formation of disks as well as guaran- 
teed their stability. 



1. Introduction 

We present results of a detailed comparison between the parameters charac- 
terizing the structural and dynamical properties of a sample of 29 simulated 
disk-like objects (DLOs) and those measured in observed spiral galaxies. 
The properties we focus on are the bulge and disk structural parameters 
and the rotation curves, as they can be constrained with available data for 
observed spiral galaxies. These data are taken principally from the com- 
pilation of galaxy structural parameters of Broeils (1992), de Jong (1996) 
and Courteau (1996, 1997). 
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Figure 1. The scale lengths and i?d (a), and the V2.2 velocities versus _Rd (b). Filled 
symbols correspond to DLOs formed in our simulations. Different symbol shapes corre- 
spond to objects formed in different realizations. Dots in (a) are from ID decompositions 
of surface brightness profiles given by Courteau, de Jong, & Broeils (1996) and Courteau 
(1997). Dots in (b) are data from the Courteau-Faber (Ha spectra) sample. Open circles 
are data from spirals with Hi rotation curves from Casertano & van Gorkom (1991). 



The DLOs have been identified in AP3M-SPH fully consistent hierar- 
chical hydrodynamical simulations, realizations of a CDM flat model, with 
f2b = 0.1, A = and b = 2.5, made using 64 3 particles in a periodic box of 
comoving side 10 Mpc (Ho = 50 km s _1 Mpc -1 ), and where an inefficient 
Schmidt law-like algorithm to model the stellar formation processes has 
been implemented (see Tissera, Lambas, k, Abadi 1997; Silk 1999). 

2. Bulge-disk decomposition 

A surface density bulge-disk decomposition was performed on the DLOs, 
using a double-exponential profile, 

E{R) = E b (0) exp[-i?/# b ] + E d (0) exp[-R/R d ]. (1) 

The resulting bulge and disk scale lengths, i?b and R^, and their ratio 
Rh/R&i are consistent with available data (Courteau, de Jong, & Broeils 
1996; de Jong 1996; Moriondo, Giovanelli, & Haynes 1999). See Fig. 1(a). 

3. Rotation curves 

The rotation curves of the simulated objects can be constructed by adding 
up in quadrature the contributions to circular rotation from the bulge and 
the disk, both formed by baryonic particles, and the dark matter halo, 

Vl(r)=Vl T {r)+VL{r)- (2) 

Our rotation curve analysis and shape modeling uses a parameterization 
based on the spatial scale R2.2 = 2.2 i?d and the velocity (mass) estimate 
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V2.2 = Kir(-R2.2)- These parameters were compared with those measured for 
the Courteau-Faber sample of bright Sb-Sc field spirals with long-slit Ha 
rotation curves, and the Broeils (1992) compilation of late-type spirals with 
extended Hi rotation curves. Deep surface photometry is available for both 
samples. In contrast to findings in other fully consistent hydrodynamical 
simulations (Navarro & Steinmetz 2000 and references therein), we find 
DLO V2.2 values that are consistent with the observational data. This is 
a consequence of disk formation with conservation of the specific angular 
momentum j. See Fig. 1(b). 

4. Luminous and dark matter contributions 

In real spirals, the relative contributions of baryons and dark matter to 
the rotation curve cannot be inferred uniquely from the observations. One 
has to postulate some hypothesis to infer the mass distribution of the halo. 
The so-called 'maximum disk (or maximal light) hypothesis' (van Albada & 
Sancisi 1986) assumes that the halo dark mass component needed to fit the 
rotation curve should be minimum, or, equivalently, that the contribution of 
the disk and bulge to the inner parts of the rotation curves should be max- 
imum. Considering this, one has (Mum/^cir)_R 2 2 — 0-85 ± 0.10, at the 95% 
confidence level (Sackett 1997). The other common constraint used is that of 
sub-maximal disks, obtained either by matching the vertical velocity disper- 
sion, scale height and scale length of a thin exponential disk which yields the 
maximum disk rotation (Bottema 1993), or, independently, showing that 
the Tully-Fisher relation of bright spirals is independent of surface bright- 
ness (Courteau & Rix 1999). Then one has (Vi mn /V cir ) R2 2 ~ 0.60 ± 0.10. 

The average relative contribution of baryons to V2.2 in our sample of 
DLOs can be measured directly. We find (Vb ar /Vd r ).R 2 2 = 0.67 ± 0.07, 
in very good agreement with Bottema (1993, 1997) and Courteau & Rix 
(1999). See Fig. 2. Most DLOs have been found to have sub-maximal 
disks. This means that the amount of baryon mass that has ended up inside 
R < R2.2 is not excessive, again as a consequence of j conservation. 

5. Conclusions 

In conclusion, the comparison between the DLOs produced in our simula- 
tions and observational data allows us to affirm that the DLOs have coun- 
terparts in today spiral galaxies. This agreement suggests that the process 
operating in Fall & Efstathiou (1980) standard model for disk formation 
(i.e., gas cooling and collapse with j conservation) is also at work in the qui- 
escent phases of the DLO formation in these simulations. However, violent 
episodes (i.e., interactions and merger events) also occur and play an impor- 
tant role in the DLO assembly (Dommguez-Tenreiro, Tissera, & Saiz 1998; 
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Figure 2. The (Vbar/Kir)fl 2 2 ratios versus the circular velocities at virial radius, V200, 
for the DLO sample. Different symbol shapes correspond to objects formed in different 
realizations. The shaded area shows the 95% confidence interval for maximum disks. 

Saiz, Tissera, &, Dommguez-Tenreiro 1999). To provide the right conditions 
for disk regeneration after the last violent episode of the DLO assembly, a 
compact central stellar bulge is needed; this will ensure the axisymmetric 
character of the gravitational potential well at scales of some kpcs at all 
times, avoiding excessive j losses in violent events. A second condition is 
necessary: the availability of gas at low z, to form the disk. The good match 
of our DLO parameters with observational data suggests that our inefficient 
star formation algorithm meets both requirements. This global agreement 
with observations also represents an important step towards making nu- 
merical approaches more widely used for the study of galaxy formation and 
evolution in a cosmological framework, i.e., from primordial fluctuations. 
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